Electronic susceptibilities in systems with anisotropic Fermi surfaces 
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The low temperature dependence of the spin and charge susceptibilities of an anisotropic electron 
system in two dimensions is analyzed. It is shown that the presence of inflection points at the Fermi 
surface leads, generically, to a TlogT dependence, and a more singular behavior, \ ~ y 3 / 4 logT, is 
also possible. Applications to quasi two-dimensional materials are discussed. 



I. INTRODUCTION 

The possible existence of quantum critical points in the 
phase diagrams of many materials has led to a detailed 
study of the low temperature behavior of the suscepti- 
bilities of electron systems. The critical properties of 
the system are determined by the energy and momen- 
tum dependence of the response function of the electron 
systeMpassociated to the order parameter in the ordered 
phasad'Btl. It has been shown that the low tempera- 
ture spin susceptibility of the isotropic electron liquid has 
an unexpected non analytic dependence on temperature, 
when high order perturbative corrections are consideredu. 
These corrections are irrelevant in the Renormalization 
Group senseffiJ'El, and do not modify the basic properties 
of the electron liquid, as described by Landau's theory. 
However, they can lead to unexpected power law depen- 
dences in many physical quantities at low temperatures, 
or change the order of the phase transitionsB. The origin 
of these non analiticitics in homogeneous response func- 
tions has been traced back to the special properties of 
2k f scattering in the isotropic electron liquids. 

It is well known that anisotropic Fermi surfaces can 
have regions where scattering becomes more singular 
than in the isotropic electron liquid, the so called " hot 
spots " . When two portions of the Fermi surface are flat 
and parallel, nesting occurs, and the susceptibilities di- 
verge logarithmically, Rex(Q,w) cx log(A/w), where Q 
is the nesting vector. A saddle point in the density of 
states leads also to logarithmic divergences in two di- 
mensions. The hot spots at the Fermi surface can be 
characterized by the frequency dependence of Imx(Q, u>), 
where Q spans the hot spots. The usual behavior in a 
Fermi liquid is Imx(Q,w) cx \u>\, in any dimension D. 
For an isotropic Fermi surface, if |Q| = 2kp , one has 
Imx(Q,w) cx \u)\( D ~ 1 ^ 2 . For D = 1 the imaginary part 
of the 2kp susceptibility approaches a constant at low fre- 
quencies. By a Kramers Kronig transformation, it can be 
shown that the real part should diverge logarithmically, 
leading to the deviations from Landau's theory which 
characterize a Luttinger liquid. 

It is also possible to show that, when Q connects 
two saddle points in an anisotropic Fermi surface, 



Imx(Q,w) 



\(D-2)/2 



This result implies the ex- 



istence of logarithmic divergences for I? = 2, which 
have been extensively studied in relation to high T c 



FIG. 1: Examples of Fermi surfaces in 2 spatial dimensions. 
Left: circular, all wavevectors of modulus 2k,F are sources 
of enhanced scattering. Right: anisotropic, wavevectors such 
as Q connecting 2 inflection points give rise to anomalous 
scattering (continuous arrow), while the rest of the FS gives 
rise to a behavior similar to the isotropic case (dashed arrow). 



superconductor and lead to deviations from Lan- 
dau's theoryllii In addition to saddle points, a generic 
anisotropic Fermi surface can show inflection points (see 
Fig.|l|). The existence of these points at the Fermi surface, 
which do not require any special fine tuning of the chemi- 
cal potential, leads tc0 Imx(Q, w) cx \ UJ \( D ~ 2 )/' 2 + 1 / i . For 
D = 2, scattering between these points is more singular 
than the 2kF scattering considered previously. 

In the present work, we analyze scattering at inflec- 
tion points in a two dimensional anisotropic Fermi sur- 
face. In the next section, we present the main features of 
the two Wop calculation, extending the method used in 
referencetl The main results are obtained in section III, 
while the finer details of the calculation are deferred to 
the appendices. Applications to Fermi surfaces of differ- 
ent shapes are given in section IV, and section V discusses 
the main results of our work. 



II. THE METHOD 

We consider a system of two-dimensional (2D) 
fermions interacting through a generic short-ranged ef- 
fective potential U (q). For the sake of simplicity, we shall 
also assume that the interaction only affects electrons of 
opposite spins, which is a reasonable approximation when 
the mciHientum dependence of U(q) is weak. It was shown 
in refsfl'tl that while the lowest order (cx U) perturbative 
corrections are well behaved, higher order corrections can 
lead to an anomalous behaviour in the low-energy prop- 
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a) 



b) 



k+Q 




FIG. 2: The second order diagrams yielding the leading tem- 
perature dependence of the susceptibility. 



erties of the Fermi liquid. To be more precise, the uni- 
form spin susceptibility of a 2D electron system shows a 
linear T dependence, which contradicts the usual Som- 
merfeld expansion in powers of (T/Ep) 2 . Such anoma- 
lous behaviour was traced back to the peculiarities of 2kp 
scattering, i.e. the occurrence of particle-hole pairs lying 
on opposite sides of the Fermi surface (FS). This special 
wavevector plays a key role in the g-dependent suscepti- 
bility of electronic systems already in the non-interacting 
case, with the appearance of a square-root singularity 
around 2kp which is directly related to the jump in the 
occupation number. If one considers the uniform sus- 
ceptibility, though, the singularities associated with 2kp 
scattering can only show up indirectly through the exci- 
tation of a virtual particle-hole pair, which explains the 
absence of anomalous corrections at lowest order in the 
interaction strength. 

In the general (non circular) case, among all the 
wavevectors connecting opposite sides of the FS, the in- 
flection points play a special role, due to the flatness 
of the Fermi surface (the extreme case being the one 
of a perfectly flat FS, or perfect nesting, which leads 
to strong instabilities). According to the previous dis- 
cussion, the second-order diagrams which lead to non- 
standard behavior are the ones containing a particle-hole 
bubble whose transferred momentum can match the spe- 
cial value Q. Such diagrams are depicted in figure |^. 
Diagrams a) and b) are vertex corrections. They have 
opposite signs and cancel in the case of a perfectly q- 
independent interaction potential: the fermion propaga- 
tors involved are the same in both diagrams, the only 
difference being in the momentum carried by the inter- 
action. To be specific, with the notations of figure |^, 
diagram a) is proportional to C/(Q) 2 , while diagram b) 
involves some momentum average of the interaction, and 
there is no reason for a perfect cancellation in the general 
case. Diagram c) is a self-energy correction, and will be 
considered separately. 



momentum limit reads: 

X(T) = 1 ,5, E/(Q) 2 A(£ P )A(£ P+Q )£(e p+Q -g p ,Q) 



(2tt) 2 (2tt) 5 



(1) 



where the Lindhard function in 2 space dimensions is 
defined as 



L(iSl,Q) = 



d 2 k n(f k ) - n(£ k+Q ) 



(2) 



(2tt) 2 *n + £k - Ck+q 

ancfi A(0 = /3/4cosh 2 (/3£/2). Such A-functions, which 
constrain momenta to lie within a shell of thickness ~ T 
from the Fermi surface, are strongly temperature depen- 
dent, and they are responsible for the leading tempera- 
ture dependence of the susceptibility. Indeed, in (|l|) we 
have omitted terms proportional to A(£ p+ Q)n(£ p ) com- 
ing from the poles of the Lindhard function, which are 
less T-dependent since they receive contributions mainly 
from regions far from the Fermi surface. Taking advan- 
tage of time- reversal symmetry (£ k = £— k)j we can write 



l(M p ,Q) 



d 2 k 



1 



1 



A£ P - A£ k A£ P 



A£ k 

(3) 



where we have defined A£ k = £ k +Q — £ k . As was pointed 
out in the previous section, the most singular contribu- 
tions to ([!]) come from regions where the momentum Q 
flowing through the Lindhard function L connects parts 
of the FS which are almost parallel, since this makes the 
denominators in eq. (^) small on large regions of fc-space. 
Otherwise stated, the scattering processes taking place 
within a particle-hole pair are enhanced around special 
wavevectors Q due to the peculiar geometry of the FS. In 
the case of a spherically symmetric FS, any momentum 
Q of modulus 2kp is a source of enhanced scattering, 
but the deviation from parallelicity is quadratic as we 
move in the direction tangent to the surface (see fig. [I], 
left). More singular is the case of inflection points occur- 
ring when the curvature of the FS vanishes, leading to 
a cubic, or even quartic dispersion (see figure |l], right), 
which is a quite generic phenomenon when dealing with 
electrons on a lattice. 

In the next subsections, we shall present the calcula- 
tion of the leading T-dependence of the diagram a) in 
the simplest circular case as well as for more complex FS 
shapes. The result for diagram b) can obtained by re- 
placing U(Q) — ► U(p — k) in eq. (fll). This can only lead 
to a change in the prefactors, but will not alter the lead- 
ing temperature dependence of the susceptibility. The 
self-energy diagram c) has a different structure, and will 
be analysed at the end of the section. 



III. UNIFORM SUSCEPTIBILITIES 



A. Isotropic Fermi surface 



After integration over Matsubara frequencies, the ver- 
tex correction a) of figure in the zero-frequency, zero- 



By choosing an appropriate coordinate system, the dis- 
persion relation around any point on a spherical FS (and, 
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generically, about non-special points of an anisotropic 
FS) can be expanded as: 



fk/« = k y +ak 2 x 



(4) 



v being the Fermi velocity at that particular point (that 
we shall identify as Q/2), and a > being related to 
the FS curvature. The above expression is assumed to 
be valid up to a momentum cutoff A which is larger than 
the one imposed by the finite temperature A-functionscil. 
Introducing q = Q Q, we can write 



6c+q/« = -(ky + q y ) + a(k x + q x y 



(5) 



were we have used the property £k+Q = £-k (reflection 
symmetry). We now change variables to 



k i — ky -\~ akj. 



(6) 



fe|| = k x + q x /2 

such that £k = vk± and 

£k+q/« = -k± -q± + 2ak\ + (3/2)ag| (7) 

Apart from the constant shift q x /2 introduced for later 
convenience, this new coordinate system is locally equiv- 
alent to polar coordinates, which would be the natural 
choice when dealing with a perfectly symmetric FS. We 
shall focus on the first term in (^|), which is independent 
of and therefore turns out to be the most singular. 
Omitting unimportant multiplicative factors, we have for 
the real part of L: 



where there is an implicit ^-function of the argument of 
the square root. We are left with a tractable expression 
for the real part of the Lindhard function, that we shall 
use to evaluate the 2-loop diagram of figure |[a. 

We now perform the remaining integrals in (Q) in the 
following order: dq±,dpj_ then dpu and dqu. The first in- 
tegral is trivial, since q± only enters in A(£ p +q). More- 
over, £ p +q is linear in q± (cf. eq. (0)) so that the integra- 
tion just gives 1/v. Thepx integral can also be performed 
straightforwardly, and we are left with an expression of 
the form 



U 2 



AU 



dq\\ J dp\\^T/v - ap 2 (11) 
d P\\ \J T / V - a P\ 



AU' 



where again we have neglected unimportant multiplica- 
tive factors and we have defined U = U(Q). Within our 
treatment, we have recovered the result that the spin sus- 
ceptibility of an isotropic 2D Fermi liquid is intrinsically 
linear in temperatureO. For a circular FS, this can be 
written as 



X{T) = Xo + XiT 



(12) 



Incidentally, our calculation suggests that the low- 
temperature correction-.to.ptbc susceptibility is positive, 
in agreement with refsflt£lliZl. 



B. Anisotropic FS with inflection points 



L = —l dk ± n(vh±) l_ A dk \\ A _ 2kx+Bk 2 
' r dk ± n(vk,) eiA - 2k±) 



v^JH 



s/A - 2k ± 



(8) 



with A = 2p± — 2ap 2 and B = 2a (the ^-function en- 
sures that the integral is real) . In the second term of eq. 
(^) we have performed the fcy integration by pushing the 
momentum cutoff to infinity. The main point is that the 
former expression can now be-integrated by parts to give 
a further A constraint on fc^El: 



1 pA/2 

L = —= I dkj_^/A- 2k A _A(vk 1 _) 
VB J -a 



(9) 



where the ellipsis stands for terms which are not confined 
to the region near the FS. By inspection of the results for 
A ^> T/v, A <§; —T/v and A w respectively, we see that 
the A(vk±) function behaves qualitatively as a 6(vk± + 
T). Therefore, to study the temperature dependence of 
the susceptibility we can replace the previous expression 
by 



L 



V 'A/2 + T/v 



(10) 



In the vicinity of an inflection point, the dispersion 
relation can be written as 



ik/v = k y - bk z x + gk 4 x 



(13) 



where b and g can be chosen to be positive. A change of 
variables similar to eq. (^|) of the previous section leads 
to 



L = 



-1 



dkj_n(vk±) / dk\ 



1 



A - 2k ± + Bkj + Ckj 



(14) 



with A = 2p^ — Sbq^p 2 — 2gpjj, B = 3bq\\ and C = 2g. 
The fcii integral can be rewritten in the form 



1 



C 1 /4(A-2fc ± ) 3 / 4 7 l + ax 2 + 



dx 



(15) 



with a = (B/yC)/y/A — 2k±. We shall be interested in 
the region close to the edge (k± ~ A/ 2), where a is large 
and positive. We can then drop the quadratic term and 
perform the integration: 



I = 



[b n (A-2ki_)Yl 2 



(16) 
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This can be integrated by parts in dk± to give 
1 f A/2 



L 



m\\ JA/2-Vql/g 



dk ± \J A - 2k±A(vk A 



(17) 



where the ellipsis stands for a term which is not confined 
close to the FS (the limits of integration account for the 
condition a > 1). Provided that the A function lies 
entirely inside the domain of integration, i.e. 



gT 

q\\ > qmin = 



the result takes the form 



y/A/2 + T/t 



(18) 



(19) 



The region of phase space we have just identified is the 
one which gives the leading temperature dependence in 
the susceptibility. Indeed, for k± outside the range of 
integration considered above (implying a < 1), the re- 
sult of the integral (|l|) is I ~ (A - 2k±_)- 3 / 4 instead of 
eq. ([l6]) , leading to a weaker (linear) temperature depen- 
dence in the final result. The same holds if we consider 
a negative q» (a < 0). 

The calculation now follows the same lines as in the 
previous case. The q± integration yields a factor 1/v, 
and the p± integration can be performed by replacing 
A(vp±) ~ S(p± — T/v)/v, which gives 



U 2 T 

'bv 4 ' 



dq\\ I dp\\ 



T/v - bq\\p 2 - gp 4 



-A 



g,l bv 4 8 V vb 2 A 2 




(20) 



Taking into account the scattering from the regions of 
the FS far from the inflection points, whose behavior is 
given by eq. (|l2|), the susceptibility reads 



x(T) = Xo + XiT-x'iTlogT 



(21) 



Once again, the sign of the correction is such that the 
susceptibility increases with temperature. However, the 
contribution coming from the other diagram b) has oppo- 
site sign. As a rule of thumb, one can argue that the over- 
all vertex-correction is positive if the effective interaction 
is peaked around Q and negative otherwise (it vanishes 
when the momentum dependence of U{q) is flat, since in 
that case the two diagrams perfectly cancel). 



C. Special inflection points 

The previous analysis assumes the existence of a 
generic inflection point along the Fermi surface. This 



is a situation which can be achieved, in an anisotropic 
system, for a finite range of values of the filling or the 
chemical potential. These points are characterized by 
the absence of a quadratic term in the expansion of the 
dispersion relation around the Fermi surface presented in 
eq. (|l^). For certain values of the parameters, however, 
which require a fine tuning of the filling or the chemical 
potential, the cubic coefficient, b, or the quartic one, g, 
in eq. ( |l3| ) can be zero as well. Two such situations are 
schematically shown in Fig|| 

We first consider the case when the cubic term in the 
dispersion relation parallel to the Fermi surface vanishes 
(6 = 0), which is realized in the t — t' Hubbard model in a 
square lattice (see point G in the top panel of figure ||), or 
in simple tight binding models on the triangular lattice, 
for instance. The susceptibility becomes more anomalous 
than in the generic case discussed previously, as can be 
seen by letting b — > in eq. (pp|). One has respectively 
A = 2p± — 3gq 2 p 2 — 2gp 4 and B = 3qq 2 . The condition 



gq 



a ^ 1 now corresponds to A — 2k± < gq'u, which modifies 
the limits of integration in eq. ( p"^ ) . Repeating the same 
arguments as before with q m i n = (T/gv) 1 ^ 4 , we obtain 



L = 



^A/2 + T/v 

v^/gq\\ 



leading to 



U 2 (T_ 
vg 



3/4 



log 



{gvA 4 ) 



(22) 



(23) 



The other possibility is that the quartic term vanishes 
(g = 0), which can occur in a tight-binding model with 
orthorombic symmetry, considering two different hopping 
parameters t a ^ tf, (see point D in the bottom panel of 
figure |). In that case, however, not only g but all the 
even coefficients in the dispersion relation vanish. This 
leads to perfect nesting between opposite branches of the 
Fermi surface, giving rise to a much more singular be- 
havior x{T) ~ logT. 



D. Self-energy correction 

After integration over Matsubara frequencies, which 
now requires some more attention due to the presence 
of two fermion lines of equal argument (a double pole in 
the complex-plane integrals), the anomalous part of the 
self-energy correction c) of figure |2| can be reduced to the 
form 



X - U 2 T / d 2 Qd 2 pd 2 k 



A(6c)A(£ P ) 
(A£ P - A£ k ) 2 



(24) 



with A£k = £k+Q — £k- The Q-integration is now re- 
stricted to the region close to (within T/v of) Q. We shall 
not go through all the calculations of the self-energy di- 
agram, which can be performed following the same lines 
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FIG. 3: Top: Fermi surfaces for different fillings in an 
anisotropic 2D system with tetragonal symmetry. Curve A 
has eight inflection points (only two are shown, Ii and I2). 
The set of all these points define a curve which goes from 
point M to point M' of the Brillouin Zone (broken curve). 
The inflection points merge in pairs when the Fermi surface 
is given by curve B. At point G in curve B, the cubic term in 
the dispersion relation parallel to the Fermi surface vanishes. 
The Fermi surface labelled C has no inflection points. 
Bottom: Fermi surfaces for different fillings in an anisotropic 
2D system with orthorhombic symmetry. Curves A and C cor- 
respond to Fermi surfaces with four inflection points. Only 
one such point in each curve is shown, Ii and I2. The set of 
all these points define a curve which goes from point M to 
point M' of the Brillouin Zone (broken curve). The quartic 
term in the dispersion relation parallel to the Fermi surface 
changes sign when going from M to M'. Thus, there is a Fermi 
surface, schematically depicted as curve B, where this quartic 
term vanishes, at point D. 

of the previous sections. The results for the temperature 
dependence are analogous to those given by equations 
(|l2|), ( pl| ) and (|2^). This can be understood by noting 
that although the denominator in eq. (p4|) is more sin- 
gular than the one of eq. (|^), the additional anomalies 
that it carries with it are cancelled by the explicit T fac- 
tor in front of eq. (p4|), leading to the same temperature 
dependence as the vertex corerction. Its sign is also the 
same as the vertex diagram a). 



The results of this section are summarized in table 



Fermi surface geometry 


X(T) 


circular 

inflection points (generic) 
special inflection (b = 0) 
nesting, saddle points 


X0 + X1T 

Xo + XiT-XiTlogT 

Xo + XiT-X3/ 4 T 3/4 logT 

Xo + Xo log T 



TABLE I: Temperature dependence of the uniform suscepti- 
bility of an anisotropic 2D Fermi liquid . The linear contribu- 
tion is always present, and is due to the portions of the Fermi 
surface away from the inflection points. The relative mag- 
nitude of the regular and anomalous contributions depends 
on the degree of flatness of the Fermi surface. The special 
case 6 = corresponds to inflection points falling on particu- 
lar symmetry lines of the Brillouin zone, and requires a fine 
tuning of the chemical potential (see text). 



IV. EXAMPLES 

A. Superconducting cuprates 

It is often assumed that a tight-binding model on a 
square lattice with nearest neighbor (t) and next-nearest 
neighbor (t 1 ) hopping reproduces well the band structure 
of the layered cuprates: 

e(k) = — 2i(cos k x + cos k y ) — At' cos k x cos k y (25) 

where t' jt w —0.25. This case corresponds roughly to the 
top set of Fermi surfaces in ng|j. The dispersion relation 
above has a saddle point at a doping by us corresponding 
to a chemical potential is Ei/hs — — The curvature 
of the Fermi surface along the diagonals becomes negative 
at a higher doping 5 C , where the chemical potential is 
E c = -8|t'| + 16^. For fillings such that E c < E F < 
Evhs, the Fermi surface has 8 inflection points. From 
these values, and the previous analysis, one can obtain 
a qualitative picture of the temperature dependence of 
the susceptibilities, when the electron density is in this 
range: 

i ) For T > \t'\, the susceptibility is determined by 
t only. As the doping is close to half filling, we expect 
X(T) cx |T|°, the result for perfect nesting. 

ii) For T < \t'\ and T > \E F - E V hs\, the same be- 
havior as in i) should be observed. 

iii) For T < \t'\ and T > \E F - E c \, the susceptibil- 
ity is dominated by contribution from the area near the 
special point discussed in section IIIC. Hence, x{T) oc 
|7f/ 4 logT. 

iv) For T < \t'\ and T < \E F - E c \, and T < \E F - 
Evhs\, the contributions from the saddle point and from 
the special point in the previous paragraph are absent. 
Thus, x(T) K 1^1 \°&(T), because of the presence of the 
inflection points. 
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FIG. 5: Correspondence between the different shapes of the 
Fermi surface discussed in the text and fillings for compounds 
of the family (TM) 2 X. 



FIG. 4: Correspondence between the different shapes of the 
Fermi surface discussed in the text and fillings for the high-T c 
superconductors. 

We can make the estimates of the crossover region 
in the T— doping plane more precise from the doping 
dependence of the coefficient of the cubic term b in 
eq.(|l^). Expanding around the saddle point, we ob- 
tain b oc \Ep — Evhs\- Hence, the crossover between 
regions ii) and iv) takes place at a temperature T* oc 
\Ep — Evhs\- Performing a similar calculation around 
the situation Ep = E c , we have b oc \f\Ep — E c \, so that 
the crossover temperature is T* oc (E c — Ep) 2 . For fill- 
ings Ep ~ E c but with no inflection points in the Fermi 
surface, we obtain a crossover to the x(T) ^ 1^1 behav- 
ior due to 2kp scattering, with a crossover temperature 
T* oc \Ep — E c \. The different regimes are schematically 
shown in Fig.||. 

Taking realistic numbers for the dispersion relation, 
our analysis predicts anomalous low temperature behav- 
ior in all the region between Svhs an d $c correspond- 
ing to the strongly overdoped region which is experimen- 
tally accessible. This shows that non-standard behaviour 
of the physical properties should be expected even in a 
regime which is usually believed to be well described by 
normal Fermi liquid theory. 

B. Quasi-ID organic compounds 

Organic conductors are often very anisotropic due to 
the planar structure of their molecules. For example, 
the salts of the family (TM) 2 X (TM=TMTTF,TMTSF 
and X=inorganic anion) are all isostructural and can be 
viewed as two-dimensional arrays of weakly coupled ID 
chains, since the electronic overlaps in the transverse di- 
rection arc 10 times smaller than in the chain direction 
(the transfer integrals in the third direction are 500 times 
smaller, and can be neglected, see for instance ref.E3). 
The band structure is well represented as: 

e(k) = — 2t a cos(k a a) — 2if, cos(fc{,&) (26) 

assuming an orthorombic structure with lattice parame- 
ters b 2a. This case corresponds to the bottom set of 



Fermi surfaces in fig. |[ The parameter t\, ~ 10 — 30meV 
sets the scale below which the FS is modulated in the 
b direction, so that the predicted enhancement of sus- 
ceptibilities due to inflection points should be observ- 
able at and below room temperature. The value of the 
anisotropy ratio r = tb /t a is large enough to ensure that 
the system is well described by a Fermi liquid picture 
down to very low temperatures. The filling factor p is 
fixed by charge transfer and varies from compound to 
compound, ranging from 1/2 to 1 hole per TM site. 
The Fermi surface has two Van Hove singularities at 
Ep = ±Evhs = ±24(1 — r), and 4 inflection points in 
all the interval < \E F \ < E VHS Taking r = 0.1, this 
corresponds to the region of fillings 0.3 < p < 1.7. In the 
absence of higher harmonics in eq.(|26|), Ep = corre- 
sponds to half- filling (p = 1), and the Fermi surface has 
perfect nesting, as e(k) = e(k + Q), where Q = (ir, tt) 
(hopping between more distant neighbors will suppress 
this effect). The points in the phase diagram where the 
topology of the Fermi surface changes, leading to differ- 
ent behaviors of the electronic susceptibility, are sketched 
in figure [|. 

In the (TM) 2 X compounds, the spin susceptibility 
shows a sizeable increase in the metallic phase up to 
room temperature (see e.g. fig. 8 of referenced) , which 
cannot be explained by "standard" Fermi liquid theory 
(the latter predicts variations on the scale of the Fermi 
temperature). On the other hand, the presence of en- 
hanced scattering close to inflection points could well be 
the underlying mechanism of this anomalous tempera- 
ture dependence, and should be taken into account when 
studying the low temperature phase transitions of such 
compounds. 

Following the same procedure used in the previous sub- 
section, the electron susceptibility will undergo a succes- 
sion of crossovers upon varying the filling, which can be 
achieved either by anion substitution or by applying pres- 
sure to the samples. 



V. CONCLUSIONS 

We have analyzed the corrections to Fermi liquid be- 
havior in anisotropic interacting electronic systems in two 
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dimensions, which arise from the existence of points in 
the Fermi surface where scattering is enhanced. Besides 
the extensively studied case of a saddle point, we have an- 
alyzed in detail the influence of inflection points, which 
do not require any special fine tuning of the chemical 
potential or the filling. The presence of these points en- 
hance the anomalous dependence on temperature which 
arise from 2k p scattering in isotropic Fermi surfacesD. We 
find that the corrections which were linear in \T\ change 
into |T|log|T|. The absence of symmetries also implies 
the lack of cancellation between different diagrams, so 
that these anomalies should be observed in both the spin 
and charge susceptibilities. 

For special fillings, more singular behavior is expected. 
In the case of systems with tetragonal or hexagonal sym- 
metry, when the Fermi surface is close to these fillings, 
the corrections to the susceptibilities go as |T| 3 ' 4 log |T|, 
showing that the existence of non integer T dependences 
does not need to violate Landau's model for the low en- 
ergy excitations of a Fermi liquid. 

We have also discussed the possible crossovers between 
the different regimes analyzed, and the experimental con- 
sequences that they may lead to. In the case of the super- 
conducting cuprates, anomalous susceptibilities should 
appear in the strongly overdoped region, above the dop- 
ing Svhs characterized by Van Hove singularities in the 
density of states. On the other hand, all of the organic 
conductors of the family (TM^X should fall in the region 
of fillings where anomalous corrections to the suscepti- 
bility are important. Of course, the analysis presented 
here should also apply to other classes of quasi two- 
dimensional systems (heavy fermion materials, Sr2Ru04, 
electrically doped 2D organic films, other organic conduc- 
tors ... ) . 

Finally, let us point out that the breakdown of the 
Sommcrfeld expansion for the spin susceptibility suggests 
that the free energy T itself has a non analytic depen- 
dence on T, once that high order interactions (2kp scat- 
texkigl-are taken into account. If, .as was proposeded 
iroll3'E£l, and numerically verified inlia, the role of tem- 
perature and magnetic field is interchangeable in the 
functional form of J-(T, H) , one can conclude that the 
anomalous T-dependences calculated here for the suscep- 
tibility are also expected in the specific-heat coefficient 
7 = C/T. 
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APPENDIX A: INFLECTION POINTS IN THE 
t t' MODEL 

We shall determine here the parameters of the disper- 
sion relation (Jl^) in the case of a tight-binding model 
on a square lattice with nearest (t) and next-nearest (t') 
neighbor hopping. Let us focus to the doping levels close 
to (5 C , the point where the inflection points of the Fermi 
surface merge in pairs on the diagonals of the BZ, lead- 
ing to the most singular corrections to the susceptibility. 
It is then natural to rewrite the dispersion relation in a 
basis rotated by 45°: 

£ = -4t(cosp cos q) + At' (cos 2 p - sin 2 q) - Ep (Al) 

where p = (k x + k y )/2 and q = (k x — k y )/2. The Fermi 
surface crosses the diagonal (q = 0) at a momentum pp 
given by Ep = —At cos pp + At' cos 2 pp. The dispersion 
relation can then be expanded as 

£ = A(p-p F ) + Bq 2 + Cq 4 (A2) 

with A = A[ts'mpp — t' s'm2pp], B — 2tcospp — At' and 
C = (4/3) [*' - (t/8) cospp]. The topology of the Fermi 
surface changes at two well-defined doping levels: 

• the curvature changes sign at a doping S — 5 C given 
by the condition B = 0. The corresponding Fermi 
energy is E c = — 8t'+16(t') 3 /t 2 and the coordinates 
of the inflection point are (p c , q c ) — (arccos2t'/i, 0), 
corresponding to the point G of figure S 

• Van Hove singularities arise at a doping 5vhs given 
by Evhs — —&t' (M-points in figure f|). 

Inflection points appear in all the region of dopings 
Svhs < S < S c , following the dashed curve of figure || 
(top panel). The consequences on the physical proper- 
ties of the system are summarized in figure ^. 

Dispersion around inflection points. The equation of 
the Fermi surface is £ = 0, which implicitly defines a 
function p = p(q). Putting the second derivative p"(q) = 
yields the coordinates (po,qo) of the inflection points. 
For S w S c , setting u = 1 — A(t' /t) 2 , we can write 

Ep — E c ( E F - i? c \ 1//2 , . . 

po= Pc + ^ N r ; ?o=(^^J (A3) 

so that the trajectory of the inflection points is parabolic 
around G. By expanding around (pcb9o)j w ^ obtain an 
equation of the form ([l3|) with 

v = Atu*' 2 ; b=-(^^-\' 2 ; g = (A4) 
t \ I2t' J ' y Atu 3 / 2 ' 
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APPENDIX B: INFLECTION POINTS IN THE 
t a - h MODEL 

We shall now derive the parameters of eq. ([l3]) for 
a tight binding model on an orthorombic lattice, with 
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anisotropic hopping (r = % /t a <C 1) . Let us rewrite for 
simplicity the dispersion relation ( |26|) as 



£ = — 2i a [cos A; + r cosp + f] 



(Bl) 



with k = k a a, p = k^b and v — Epjlt. The equa- 
tion of the Fermi surface is k = arccos(— v — rcosp). 
The number of electrons per site is approximately given 
by p — 2tt~ 1 arccos(— v). The Fermi surface has 2 Van 
Hove singularities at Eyas — ±2i(l — r), and 4 inflec- 
tion points for any < \Ep\ < Eyus- At half filling 
{Ep = 0), the two branches of the open Fermi surface 
are perfectly nested. The physical consequences of the 
changes in the Fermi surface topology occurring at those 
special fillings are sketched in figure 0. 

Dispersion around inflection points. By setting 
k"(p) = we find that the inflection points are located 
at 



Po 



arccos 



; fco = arccos [ — v — 



1 



(B2) 



For filling levels close to v = (half filling), the location 
of the inflection points describes a straight line of slope 
— t j\J\ — v 2 in the (k,p) plane (see bottom panel of fig- 
ure After a rotation of the coordinate axes, we obtain 



an equation of the form (13) with 



2i\A - v 2 ■ b = 



6VT 



; g 



24(1 - V 2 fl 2 

(B3) 



The approach to the perfect nesting situation at half fill- 
ing is signalled by a vanishing g, the coefficient of the 
quartic term in the dispersion relation. 
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